import numpy as np
import matplotlib.pylab as plt

d = np.loadtxt('solar_panel_ang_25_45.dat')

plt.figure(figsize=(10,8))
ax = plt.subplot(111,projection='mollweide')
ax.grid(True)

x = []
y = []
for i in range(100000):
    x.append(np.radians(d[i,1]))
    if d[i,2]>180:
        y.append(np.radians(d[i,2]-360))
    else:
        y.append(np.radians(d[i,2]))

ax.scatter(y,x,s=0.01)
plt.show()

